import numpy as np

sample_list = [50,60,75,85,100,30,48]
sample_num = len(sample_list)
x = np.power(np.prod(sample_list),1/sample_num)
print(x)

#方法2
from scipy.stats import gmean
y=gmean(sample_list)
print(y)

#方法3
def prod_f(x,y):
    return x*y

from functools import reduce
z = np.power(reduce(prod_f,sample_list),1/sample_num)
print(z)